Arctangent Integral

Function: double gsl_sf_atanint (double x)
Function: int gsl_sf_atanint_e (double x, gsl_sf_result * result)
These routines compute the Arctangent integral AtanInt(x) = \int_0^x dt \arctan(t)/t.